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The momentum or velocity autocorrelation function C{t) for a tagged oscillator in a finite har¬ 
monic system decays like that of an infinite system for short times, but exhibits erratic behavior 
at longer time scales. We introduce the autocorrelation function of the long-time noisy tail of C{t) 

(“a correlation of the correlation”), which characterizes the distribution of recurrence times. Re¬ 
markably, for harmonic systems with same-mass particles this secondary correlation may coincide 
with the primary correlation C{t) (when both functions are normalized) either exactly, or over a 
significant initial time interval. When the tagged particle is heavier than the rest, the equality 
does not hold, correlations shows non-random long-time scale pattern, and higher order correlations 
converge to the lowest normal mode. 

PACS numbers: 05.40.Ca, 05.20.Gg 


I. INTRODUCTION 

The theme of fluctuations in finite systems of harmonic 
oscillators emerges naturally in both application and the¬ 
ory. From a theoretical point of view, the study of the 
stochastic dynamics of a tagged degree of freedom in fi¬ 
nite harmonic systems provides a valuable illustration, 
and often more than that, of the role of the thermody¬ 
namic and weak-coupling limits, ergodicity, thermaliza- 
tion, recurrences, synchronization, and other basic con¬ 
cepts in nonequilibrium phenomena [IH3- Another rele¬ 
vant area is Langevin dynamics generated by a coupling 
to a finite harmonic bath(s), and its application to meso¬ 
scopic systems and networks; see [9l-[lq. 

Being nonergodic, the capability of harmonic systems 
to illustrate general phenomena in statistical mechanics 
might seem doubtful at first glance. By means of a canon¬ 
ical transformation a harmonic system of any size can be 
transformed into a collection of independent oscillators, 
or normal modes; since the energies of normal modes are 
the integrals of motion, a single isolated harmonic system 
does not equilibrate and is not very interesting from the 
point of view of statistical mechanics. 

A more fruitful approach is to consider an ensemble of 
harmonic systems, assuming that in the past they were 
in contact with a larger thermal bath in equilibrium at a 
given temperature, and that the initial normal modes of 
the ensemble are distributed according to the canonical 
distribution. Within this framework, one evaluates sta¬ 
tistical averages of dynamical variables over the ensem¬ 
ble of the system’s initial coordinates rather than over 
time. Such averages show the transition of the ensem¬ 
ble to thermal equilibrium in the limit of a large num¬ 
ber of particles, and thus the nonergodic nature of har¬ 
monic systems does not explicitly manifest itself, and for 
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most cases is inessential. It should however be stressed 
that this framework, which is standard for most works 
on stochastic dynamics of harmonic systems both clas¬ 
sical and quantum, assumes a very special type of cou¬ 
pling between the system and the external thermal bath: 
this coupling justifies the initial conditions for the the 
system’s degrees of freedom, yet is assumed to be suffi¬ 
ciently weak, or completely turned off, as not to affect 
the system’s further dynamics. 

The inequivalence for harmonic systems of ensemble- 
and time-averages, together with the almost exclusive ex¬ 
ploitation in literature of the former, does not necessarily 
entail that the latter are inadequate. Rather, we intro¬ 
duce in this paper a new class of time-average correlations 
(we call these secondary correlations) which characterize 
recurrences in finite harmonic systems. For systems of 
same-mass particles, these correlations are shown to be 
very close, and under certain conditions exactly identical, 
to the conventional (primary) time correlations defined 
with ensemble averaging. This implies that for finite non¬ 
ergodic systems, the use of both ensemble and time av¬ 
erages may give meaningful complementary descriptions, 
and that correlations with the two types of averaging may 
be related in some subtle way. 


II. SECONDARY CORRELATIONS 

Consider the temporal autocorrelation function 
(A(O)A(t)) of a dynamical variable A in a finite system of 
size L - typically, such a function exhibits two distinctive 
regimes, separated by a crossover time tc of order L/v, 
where v is the speed of signal propagation in the system. 
For short times t < tc, the variable does not feel the 
presence of the boundaries, and the correlation function 
decays in a smooth regular way, following the same laws 
as for an infinitely large system. On the other hand, for 
longer times t > tc the dynamics of the variable are af¬ 
fected by signals reflected from the boundaries. For long 
time regimes such as this, rather than decaying smoothly 



2 


the correlation functions may exhibit erratic, apparently 
noisy, behavior 0,i. 

We illustrate this behavior in Fig. [T] by way of the nor¬ 
malized momentum correlation function for the central 
particle in a harmonic chain with fixed ends. The Hamil¬ 
tonian of the system is 
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which describes iV -|- 2 linearly coupled particles indexed 
1 = 0,1,..., -|- 1 with terminal particles fixed, with dis¬ 

placement go = Qn+i = 0. Assuming N is odd, the 
middle particle indexed 
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has normalized (C'(O) = 1) momentum correlation func¬ 
tion 
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where the prime indicates that the summation is only 
over odd j. In this expression (we outline its derivation 
in the Appendix), the ujj terms are frequencies of normal 
modes 


ujj = 2(jj sin 


2{N+1) 


( 4 ) 


where w is the frequency of a single oscillator, and the 
average (• • •) is taken over the equilibrium ensemble of 
initial conditions. For t < tc the correlation function 
Ci(t) is very close to that of an infinite chain, given by 
the Bessel function 
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FIG. 1: The normalized momentum correlation function Ci{t) 
for the central particle (i = io = 51) of the harmonic chain 
with fixed ends with Hamiltonian o with N = 101, given by 
Eq. ®. The time unit is l/2tJ. The time of the crossover 
from the regular dissipation to “stochastic” regimes is tc ~ 
200 . 


In this paper we propose to characterize the irregular 
part of the function Ci(t) in another way, which is more 
in the spirit of nonequilibrium statistical mechanics than 
the mathematics of almost periodic functions. Namely, 
observing that for large t the correlation function Ci{t) 
appears to behave like stationary noise, we are encour¬ 
aged to characterize it by a new correlation function 

A(i) = H 

Wi ('r))r 


Ci{t) Ri Coo(i) = lim Ci{t) = Jo{2ujt). (5) 

i,N—^oo 

This can be readily justified by approximating the sum 
([3]) with an integral, and recognizing the latter as the 
well-known integral representation of Jo(2a;t), see e.g. 

[II- 

More interesting from the perspective of this paper is 
the regime t > tc in which the correlation Ci{t) becomes 
irregular, see Fig. 1. It can be shown that the function 
Ci{t) given by ([3]) belongs to the class of almost periodic 
functions: any value c which the function achieves once 
is achieved again, infinitely many times. Traditionally, 
such functions are characterized by the average frequency 
with which they return to c, or by the reciprocal, i.e. the 
mean recurrence time t(c). For correlations of type ® 
with large iV, the famous result for the recurrence time, 
first obtained by Kac Q (see also Qai), 

r(c) ~ (6) 

implies that recurrences of order c ~ or larger are 
exponentially rare. This result resolves, or rather (be¬ 
ing derived for a model system) shows the direction of 
resolution for the paradoxes of irreversibility [I|. 


defined with the time average 


(-)t = lim ^ [ {...)dT. 
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Since we are only interested in the interval t > tc when 
Ci{t) behaves irregularly, one might prefer to set the 
lower integration limit in definition ([5]) to tc instead of 
zero. However this would only be an unnecessary com¬ 
plication, as the limit T —> oo makes the two definitions 
numerically equivalent (assuming always that the inte¬ 
gral from 0 to tc converges). 

We shall refer to Di{t), defined by relations ([7]) and 
® , as the secondary correlation function, and call Ci (t) 
the primary one. We would like to promote the secondary 
correlation Di{t) as a meaningful statistical tool for char¬ 
acterizing the distribution of recurrences times in a sys¬ 
tem of finite size. Such information is not contained in 
the Kac formula (l6|) for the average recurrence time r, so 
the two functions r(c) and Di(t) do not duplicate each 
other but describe recurrences in complementary ways. 
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III. RELATION TO PRIMARY CORRELATIONS 


Since the primary and secondary correlations Ci{t) and 
Di{t) characterize recurrences at different levels and are 
defined using different types of averaging (over ensemble 
and time, respectively), the existence of any specific rela¬ 
tion between them is perhaps a priori unexpected. Yet a 
simple numerical experiment with Eqs. (1-6) suggests, for 
the middle atom of a chain with fixed ends, the equality 

C.„(t)=Ao(t). (9) 


Closer scrutiny reveals that the equality is exact and 
holds for any t, such that the secondary correlation com¬ 
pletely repeats the structure of the primary one for both 
regular [t < tc) and noisy {t > tc) domains and has the 
same crossover time The proof follows immediately 
from the relation 


{cOSOJjT COSUJjiir + t))r 


Sjf 

= —— cos Wot 
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which holds for an arbitrary spectrum of (nonzero) nor¬ 
mal mode frequencies {ujj} and can be verified by direct 
evaluation (with the help of L’Hospital’s rule). For t = 0 
this may further be reduced to the familiar orthogonality 
relation for the Fourier basis, and thus can be considered 
a generalized form of the latter. From (l3|) and m one 
obtains for the non-normalized secondary correlation 


(Cioir) C,„{t + t))r = 
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Normalizing this function to unity at t = 0 by dividing 
it by 


{Cl{T))r 


If 2 y N + 1 
2 \N+l) 2 


1 
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one obtains the normalized secondary correlation 
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FIG. 2: The primary momentum correlation function Ci{t) 
given by Eq. (I14II (solid line) and secondary correlation func¬ 
tion Di{t) given by Eq. (I17II (dashed line) for the harmonic 
chain with Hamiltonian © with N — 101 for particles i = 20 
(top), i = 30 (middle), and i = 40 (bottom). The difference 
between Ci{t) and Di{t) becomes noticeable for t > to where 
to depends on i non-monotonically: to ~ 80,120, 90, from top 
to the bottom. 


function for particle with arbitrary index i has the form 
(see Appendix) 

N j —2— ■ 

C,{t)=Y,Alcosujjt, Ay = ^^^^sin-^^. (14) 

For the middle particle i = io = {N + l)/2 this is reduced 
to (O, whereas for the other particles normal modes en¬ 
ter the expression (03) with different amplitudes . As 
one can immediately verify, the exact equality of primary 
and secondary correlations does not hold in these cases. 
An important example when this equality does hold for 
any particle is a harmonic chain with periodic bound¬ 
ary conditions. In this case the momentum correlation 
for each particle is a superposition of equally weighted 
normal modes na 

1 ['^j\ 

= A ^iv ) ’ 
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cosujjt (13) 


which coincides with the primary correlation Ap (t), 
Fq.®. 

One may observe that for the above derivation it is es¬ 
sential that the primary correlation Ci{t) takes the form 
of a superposition of cosines with equal weights, as in 
Fq. ®. In general this, of course, is not the case. For 
example, for a chain with fixed ends described by the 
Hamiltonian o. the normalized momentum correlation 


and repetition of the above derivation leads again to the 
exact equality Ci{t) = Di{t) for any particle of the sys¬ 
tem. 

So far, even with the above examples of its validity, 
the equality of primary and secondary correlations may 
appear as no more than a curious coincidence. How¬ 
ever, further numerical exercises reveal that even when 
equality does not hold exactly, it remains a very good 
approximation for the initial time interval t < to, see 
Fig. 2. The duration of this interval, to, is found to 
depend non-monotonically on particle position i, and for 
any i be equal or shorter than the crossover time, to < tc- 
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Respectively, for t < tg both primary and secondary cor¬ 
relations coincide with the primary correlation for the 
infinite chain, 


Di{t) = C,{t) = Caoif) = '/o(2wi), t <tg< tc- (16) 

The proof of the approximate equality (1161) can be car¬ 
ried out as follows. From the expression (fHl) for Ci{t) 
and the definition (l7|) for Di, and using the relation (TTOl) . 
one gets 


D^{t) = cosujjt, (17) 

E 4 

i=i 

or taking into account the expression (|3]) for normal mode 
frequencies 
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Recognizing here the generating function for Bessel func¬ 
tions 


cos(a;sin0) = Jo{x) +2''^J2k{x) cos(2fc0), (19) 

k=l 

Di(t) can be written as a superposition of Bessel func¬ 
tions 


Di{t) — Jo{2ijjt)+ ''^^ Sik (20) 

k=l 

with coefficients 

s«=v^i:4™(j^). (21) 

E 4 

A simple analysis of this expression shows that given i, 
the coefficients Sik are nonzero only for five sets of k: 
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where s = 1, 2, 3,... Note that this expression is invari¬ 
ant under the transformations i —>■ {N -|- 1) — i, reflecting 
the symmetry of the left and right sides of the chain. One 
can observe that for large N and i not too close to the 
end or to the middle of the chain the coefficients Sik are 
nonzero only for large indices k. For instance, for the 
chain with N = 101 and the particle i = 20, coefficients 
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FIG. 3: The secondary correlation functions D 3 o(t) (top) 
and (bottom) according to the exact expression (1171) 

(solid lines) and the approximation (1231) (dashed lines). The 
insets show apparently random behavior of Dzo{t) and D^oit) 
at longer times. 


Sik are nonzero only for fc = 40,80,124,... As a result, 
for t not too large in the expression (|20)) , the dominating 
contribution comes from the first term Jg(2Lot)^ while the 
corrections given by the sum E^i J 2 k{‘^xit) involve 
Bessel functions of large orders which are negligibly small 
for a significant time interval t < tg [l^ . 

The above consideration not only justifies the equality 
Di{t) = Ci{t) = Jg{2ujt) for t < tg, but also accounts for 
a curious non-monotonic dependence of tg on the tagged 
particle index i, which we noticed empirically in Fig. 2. 
For example, according to (1^^ . for N = 101 and par¬ 
ticles i = 20,30,40 the minimal indices k for which Sik 
takes nonzero values (—4/3, —4/3,1/3) are k = 40, 60,44, 
respectively. Then keeping only the leading and first cor¬ 
rection terms in the exact expression (EOD, one gets 

D2o{t) = •A)(2wt) — - J8o(2wt), 

Dsoit) = •/o(2wt) — - Ji2o(2(jJt), 

Dioit) = •A)(2wt) + — J88(2wt). (23) 

One can verify that these approximations describe the 
initial deviation of Di{t) from Caa{t) = Jg{2ujt) very well 
indeed (Fig. 3 shows this for particle i = 30 and i = 40). 
Since for small arguments Ji{x) decreases with order i, 
it is clear from (051) that the second correction terms for 
particles i = 20,40 involve Bessel functions of smaller 
orders, and thus become essential at earlier times than 
for particle * = 30. 

If one applies a similar analysis to the primary corre¬ 
lations Ci{t), Eq. (fT4l) . one gets a familiar approximate 
relation for the left side of the chain [^, |^, [l^ 

C'i(t) = Jg{2,ujt) — Jii{2ujt). 


(24) 
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Here, in contrast to corresponding relations (|23)) for 
Di{t), the order of the second Bessel function, which de¬ 
scribes effects of finite size, increases monotonically (lin¬ 
early) with particle index i, and so does the crossover 
time tc- 

In order to study the dependence of the characteris¬ 
tic time to I during which Ci{t) = on i in a more 

quantitative way, let us consider the function 

6,it) = Qit)-D,it), (25) 

which is zero when the two correlations coincide for t < 
to and fluctuates at longer times. For a given t, let us 
define to somewhat arbitrarily as the time at which 6i{t) 
reaches its first local minimum or maximum; see Fig. 
4(a). Similarly, we can define the crossover time tg as 
the moment when the function 


A*(t) =a(t)-Coo(t) (26) 


has its first local extremum, recalling that Coo(t) = 
Jo(2w) is the correlation in an infinite system. Using 
these definitions, we record observations of to and for 
N = 101 in Fig. 4(b), as a function of particle index i. 
Whereas tc increases linearly as we approach the central 
particle, to coincides with tc for t < U = 34 and linearly 
decreases for i > ii. As we already know, the primary 
and secondary correlations coincide for the middle parti¬ 
cle to = 51, so Sig{t) is identically 0 and to diverges here. 
Somewhat unexpectedly, we find that to also diverges, 
i.e. C'i(t) = Di{t) identically, for ii = 34 (and of course 
the symmetric case i 2 = {N -|- 1) — U = 68). Therefore 
it would appear that to{i) diverges whenever it changes 
from increasing to decreasing, or vice versa. Further cal¬ 
culations for different N show that in general the exact 
equality Ci{t) = Di{t) holds for particles with indices 


*0 = 


N + 1 


*1 = 


Af-bl 


*2 = 


2{N + 1) 


(27) 


provided of course that these expressions are integers. 
For N = 101 there are three such particles (*o = 51, ii = 
34,12 = 68), two for N = 200 (U = 67,12 = 134), and 
none for N = 100. 

Let us show that this phenomenon is readily ac¬ 
counted for with Eas. (l20l) - (l2^ for the secondary correla¬ 
tion Di{t). First, from inspecting (I22|) one might observe 
that for i < io the minimum k for which Sik is non-zero 
is k = 2i and comes from the set k = 2{N + l)(s— 1) + 2i 
with s = 1. This yields the approximation 


Di{t) = Jo{2ujt) - ^ J4j(2a;t), (28) 

which we already used for D 2 o{t) and Il 3 o(t) in (051) . 
It differs from the corresponding approximation (1241) for 
Ci{t) only by the factor 4/3 in the second term. Then, 
from (05|) and (051) , the difference functions defined above 
by relations (05l) and (051) take the form 


^i{t) = 3 •hi{2u}t), 


Ai(t) = -Jii{2ujt). 


(29) 
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FIG. 4: Top plot (a): The functions Ai{t) = Ci(t) — Jo(2a;t) 
(solid line) and Si{t) = Ci{t) — Di{t) (dashed line) for particle 
i = 40. The characteristic times U and to are defined as 
times at which Ai{t) and Si{t), respectively, have their first 
local extremum. Bottom plot (b): The characteristic times U 
(x) and to (-I-) for particles with indices i < 51 for the left side 
of the chain described by Hamiltonian ([T]) with N — 101. For 
particles io = 51 and U = 34 the time to diverges (illustrated 
by an arrow pointing upward), indicating the exact equality 
Ci{t) = Di{t). 


Since these two functions have local extrema at the same 
time, by definition we have to = U- Furthermore, since 
the position of the first maximum of the Bessel function 
Ji{t) increases approximately linear with i [ 13 : Eq. (1551) 
explains the equality of the characteristic times tdi) = 
to^i) and their linear increase for f < ig in Fig. 4(b). 

As i gets larger still, one observes from (l22l) that a 
minimal k for which Sik yt 0 is /c = 2{N -|- 1) — 4z and 
comes from the set k = 2{N + l)s — 4i with s = 1. 
In this case for Di{t), instead of (1551) . we have another 
approximation 

Di{t) = Jo{2u!t) + — Jci{2ujt), 

a = AiN + l)-8i, (30) 

which we already used for D 4 o{t) in (1551) . Since the pri¬ 
mary correlation Ci{t) is still given by (1551) . the difference 
function 5i{t) = Ci{t) — Di{t) in this case reads 

Si{t) = Ja(2u}t) - Jn(2wt) ~ Jc,i2wt). (31) 

The position of its first extremum increases approxi¬ 
mately linearly with a [13 and, as follows from dsni, 
linearly decreases with i. This explains the behavior of 
to(*) for i > U in Fig. 4(b). 

The transition of to(0 from a positive to a negative 
slope (over the domain [0, 51]) occurs at * = i\, for which 
^2 = 2{N + 1) — Ai (the minimal value of the set k = 
2{N -b l)s — 4i) becomes less than or equal to ki = 2i 
(the minimal value of the set k = 2{N + l)(s — 1) -|- 2i). 
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Then the equality fci = ^2 gives ii = {N + l)/3, which is 
consistent with our empirical findings (E71) . 

The exact equality Ci{t) = Di{t) for i given by (l27l) 
can be readily verified using the following expression for 
the primary correlations 

OO 

Ci{t) = Jo(2uit) + Tik J2fc(2wt), (32) 

with coefficients 

life = 2 ^ Afj cos ■ (33) 

i=i ^ ^ 

These relations are similar to (1201) and (1211) for Di{t) and 
can be derived in a similar way 0- For i = io,*i ,*2 
given by dUD, one can verify directly from (1331) and m 
that Sik = Tik for any k. Then the comparison of (1321) 
and (1201) gives for those values of i the exact equality 
C,{t) = D,{t). 


IV. HEAVY IMPURITY PROBLEM 
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FIG. 5: The primary (solid line) and secondary (dashed 
line) momentum correlation functions, C{t) and D(t), for the 
heavy impurity problem described by Hamiltonian dMl) with 
mass ratio p = m/M = 0.1 and N = 50. The top, middle, 
and bottom figures show the evolution of the correlations on 
short (t < tc), intermediate (t ~ tc), and long (t tc) time 
scales. The inset shows correlations on the short time scale 
for an impurity that is twice as heavy; p = 0.05. Time is in 
units of l/2a;. 


So far we have discussed finite harmonic systems of 
similar particles. If a tagged particle is heavier than 
the rest, it turns out that the equality of primary and 
secondary correlations, C{t) and D(t), does not hold. 
Though structurally similar - D{t) looks like a coarse¬ 
grained copy of C{t) - the two correlations are quite dis¬ 
tinctive on any time scale; see Fig. 5. In particular, 
the approximation of exponential relaxation for t < tc, 
while good for C'(t), is noticeably worse for D{t). An¬ 
other observation is that for t ^ tc both correlations, 
being apparently random on a short time scale, show on 
a larger scale a noisy yet periodically repeating pattern; 
see the bottom plot in Fig. 5. This feature, absent in 
systems of equal-mass particles, is made all the more ob¬ 
vious when considering higher order correlation functions 
Ck{t), defined recursively as 


Ck+i{t) 


{Ck{T)Ck{T + t))r 

{Cl{r))r 


(34) 


assuming new notations for C{t) = Ci{t) and D{t) = 
C 2 {t). (In this section we use the notation Ck{t) with a 
subscript referring to the correlation order, rather than 
to the index of a particle). For the heavy impurity prob¬ 
lem one finds that as the order k increases the apparent 
randomness of correlations Ck (t) on the time scale t > Tc 
quickly diminishes^ and higher correlations converge to 
the normal mode with the lowest eigenfrequency Hi: 


Ck{t) —>■ cos(Uit); (35) 


see Fig. 6. Below we outline a theoretical framework 
underlying these empirical observations. 


Consider a cyclic chain of 2N particles of mass m and 
an impurity of mass M > m described by the Hamilto¬ 
nian 


H = 


2M 


2N 2 

A. 

2m 


E 


mu) 


2N-1 


E “ ®+i)" 


2 = 1 


muj 


liQ - 9l)^ + (Q - 92Ar)^], 


(36) 


where P and Q are the momentum and coordinates of 
the impurity. Using a diagonalization method similar 
to that described in the Appendix (see Q for details), 
one can show that the normalized momentum correlation 
function for the impurity C{t) = {P{0)P{t))/{P‘^{0)) is 
again an almost periodic function, now of the form 


2Af-l 


C{t) = Ci{t) = ^ Aj cosfljt. (37) 

The amplitudes Aj in this expression are given by 


2Af-l 


..S- 


(38) 


where 


u!j = 2ijj sin 


ei = —2^2 


2{2N+1)’ 
2 


2N + 1 


sm ■ 


2V + 1’ 


(39) 
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time 


FIG. 6: The primary correlation Ci{t) (solid line), and two 
higher correlations Czit) (dashed line) and C^it) (dotted line), 
defined by Eq. (OH) . at long time t^tc for the heavy impu¬ 
rity problem with fi = 0.1 and N = 50. Higher correlations 
converge to the lowest normal mode cos Hit (dot dashed line). 

and /t = m/M is the mass ratio. 

Due to the system’s symmetry only the modes with 
zero and odd indices contribute to the superposition (1371) . 
Their frequencies flj (j = 0,1, 3, • • • 2N — 1) for M ^ m 
cannot be expressed in closed form and must be evaluated 
as roots of the secular equation Q 

2N-1 2 

G{z) = — 2 — — - — 5 -= 0. (40) 

This transcendental equation has N + 1 solutions z = 
ilj, j = 0,1,3, •••2iV — 1. It can be verified that one 
solution is the zero frequency Dq = 0 , which reflects the 
translational invariance of the system. The remaining N 
nonzero roots Di, D 3 ,... n 2 N-i he in the interval ( 0 , 2a;) 
and must be evaluated numerically. 

With the set of eigenfrequencies Dj found, one may 
calculate the amplitudes Aj with (l38l) and evaluate the 
primary correlation C(t) by carrying out the summation 
in dST]). Then, using m, for the secondary correlation 
D{t) one obtains 

f 1 1 

D{t) = C2it) = C2 IaI + - AjCosOjtj. (41) 

with normalization coefficient 

Fig. 5 presents C{t) and D{t), calculated with Eqs. (1571) 
and m, for /i = 0.1 and N = 50. In a similar man¬ 
ner, one can obtain the expression for order-fc correlations 


from (l34l) 

f / A \°"‘ 1 

Ck{t) = Cfc < ^ 0 " + cosfljt \ (43) 

with powers ak = 2 ^“^ and normalization coefficient 



While expression (H51) is a superposition of -I-1 modes, 
one can observe that for larger k the main contribution 
comes from the mode with eigenfrequency Oi, such that 
Ck{t) ~ cos Hit. This can be accounted for by notic¬ 
ing that the sequence of coefficients {Aj} has Ai as its 
maximum element and is monotonically decreasing for 
j > 0 . For example, for ^ = 0.1 and = 50 we find that 
^0 = 0 . 09 , Ai = 0 . 17 , A 2 = 0 . 15 , ^3 = 0 . 12 , ... (approx¬ 
imately). For the primary and secondary correlations 
involving Aj and A^ respectively, such an insignificant 
difference in values hardly plays a role. But for higher- 
order correlations the maximum of the set {} (still at 
j = 1 ) may be orders of magnitude greater than any other 
element. As a result, the superposition in (H5)) is increas¬ 
ingly dominated by the term with , and higher-order 
correlations quickly converge to the first normal mode; 
Cfc(t) ~ cos Hit. 

For systems of same-mass particles the set of normal 
mode amplitudes, given by the second equation in (I14L 
is a periodic function of the mode index j and has no 
single maximum. In this case the reduction of higher- 
order correlations to a dominating normal mode does not 
occur. 


V. CONCLUSION 

Temporal autocorrelation functions (A(O)A(t)) are of¬ 
ten evaluated in the thermodynamic limit, in which case 
they typically decrease in a regular (non-random) fash¬ 
ion, either monotonically or non-monotonically. In fi¬ 
nite systems, autocorrelation functions themselves be¬ 
come noisy at long time scales t > tc\ this illustrates 
recurrences in the dynamics of the tagged variable due 
to reflections of sound off boundaries. In this paper we 
introduced and studied some properties of the secondary 
correlation function D{t) defined as an autocorrelation 
function of the primary correlation function C{t). If it ex¬ 
ists, the characteristic time of decay for D{t) determines 
the time-scale of a typical “period” for C{t), which in 
turn may be associated with the typical recurrence time 
of the targeted variable. These “typical” times may how¬ 
ever be ill-defined mathematically (as is indeed the case 
for the harmonic systems discussed above), so to be more 
precise the secondary correlation D{t) can be described 
as a function characterizing a distribution of recurrence 










times: for a given t, a larger value for D(t) corresponds 
to a greater probability (density) that C{t) will return to 
an assigned value in time about t. Comparing the sec¬ 
ondary correlation D(t) with the mean recurrence time 
t(c), Eq. ([6|), the latter being more prevalent in litera¬ 
ture, one notices that the two functions give complemen¬ 
tary descriptions: while r(c) characterizes the number of 
returns to an assigned value c, the secondary correlation 
D{t) gives the distribution of return times regardless of 
the assigned value of c. 

One interesting result is the equality C{t) = D{t) 
for systems of same-mass particles. The equality is ei¬ 
ther exact for all t or a very good approximation over 
the initial interval t < to whose duration to depends 
on the tagged particle’s position non-monotonically. Al¬ 
though its derivation is quite simple, the equality of pri¬ 
mary and secondary correlations may be a remarkable 
property, especially considering that the former is de¬ 
fined over the ensemble and the latter with time av¬ 
eraging. We restricted the discussion to the simplest 
case of one-dimensional harmonic systems, but an exten¬ 
sion to higher dimensions appears to be straightforward. 
Whether the equality, or perhaps some other relation, be¬ 
tween primary and secondary correlations still holds for 
nonlinear systems is an open question. 

Like the primary correlation, for long time-scales the 
secondary correlations also develop noisy tails (see the in¬ 
sets in Fig. 3), which themselves can be characterized by 
correlations of higher order. In turn, this new tertiary 
function has the same structure as the secondary and 
primary functions, exhibiting regular decay over shorter 
times and fluctuating over longer times. Thus one can 
construct an infinite hierarchy of higher order correla¬ 
tions whose order-scaling properties are also interest¬ 
ing to study. Of course, for the particular cases when 
the equality D{t) = C{t) holds exactly for all t, e.g. a 
harmonic chain with periodic boundary conditions, all 
higher order correlations are identical. We have studied 
higher correlations in the context of the heavy impurity 
problem. Here the equality of primary and secondary 
correlations does not hold, and the primarily correlation 
displays a non-random idiomatic pattern on long time 
scales, which becomes even more visible in correlations 
of higher orders. Indeed the sequence of higher-order 
correlations converges to the lowest normal mode. 
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Appendix 

In this appendix we outline the derivation of expres¬ 
sions and m for the momentum correlation func¬ 


tions of i-th particle in a harmonic chain with fixed ends, 
described by the Hamiltonian ©• Using the normal 
mode transformation 

^ AT N 

9* ~ ~/m Pi ~ ^ijPj 

i=i 

with coefficients Aij given by m and taking into ac¬ 
count the orthogonality relation the 

Hamiltonian © is diagonalized into the form of uncou¬ 
pled normal modes 

1 ^ 
i=i 

with frequencies ujj given by ©• The normal modes are 
governed by the Hamiltonian equations 

dH „ . dH 

and evolve as 

Pj (t) = Pj (0) cos (jjjt — ujjQj (0) sin ojjt, 

Qj{t) = Qj{0)cosujjt + u}~^Pj{0)smu}jt. 

Assuming that initially the system is in equilibrium 
with canonical distribution function pe = , 

correlations of the normal modes’ initial values are 
(P,(0)P,.(0)) = and (Q,(0)p,.(0)) = 0. Then 

s 

= (Py(O)Pj(O)) COSujjt = cosujjt 
and the momentum correlation of the i-th particle is 

N 

{p^{0)p^{t)) = m ^ AyAyv(P^-,(0)Pj(t)) 
id'=i 

N 

m .9 

= ^ Z^^^jCOSW.P 

Division of this expression by (p^(0)) = m/fi gives the 
normalized correlation function (imi . In the case of the 
middle particle i = {N + l)/2 (assuming N is odd), = 
2/ {N -I- 1) for odd j and zero otherwise. In this case, one 
obtains the normalized correlation function in the form 
The correlation (1151) corresponding to the periodic 
boundary condition can be derived in a similar way. For 
the extension to the heavy impurity problem see e.g. Q. 
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